Blueshift of the surface plasmon resonance in silver nanoparticles studied with EELS 
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We study the surface plasmon (SP) resonance energy of isolated spherical Ag nanoparticles in 
the diameter range 3.5-26 nm with monochromated electron energy-loss spectroscopy. A significant 
blueshift of the SP resonance energy of 0.5 eV is measured when the particle size decreases from 26 
down to 3.5 nm. We interpret the observed blueshift using three models for a metallic sphere em- 
bedded in homogeneous background material: a classical Drude model with a homogeneous electron 
density profile in the metal, a semiclassical model corrected for an inhomogeneous electron density 
associated with quantum confinement, and a semiclassical nonlocal hydrodynamic description of 
the electron density. We find that the latter two models provide a qualitative explanation for the 
observed blueshift, but the theoretical predictions show smaller blueshifts than observed experimen- 
tally. We attribute the larger measured blueshift to effects due to the presence of the substrate, 
which breaks the symmetry of a particle embedded in a homogeneous medium. 



INTRODUCTION 

Surface plasmons are collective excitations of the elec- 
tron gas in metallic structures at the metal/ dielectric 
interface The ability to concentrate light with 

SPs 0] and to enhance light-matter interaction on a sub- 
wavelength scale enables few- and even single-molecule 
spectroscopy when the size of the metallic structures 
is decreased to a few nanometer These collective 
excitations are usually well-described by the classical 
Drude model for nanoparticles with dimensions of tens 
of nanometer and larger In the quasistatic limit, 
i.e. when the wavelength of the exciting electromagnetic 
wave considerably exceeds the dimensions of the struc- 
ture, the local-response Drude theory predicts that the 
resonance energy of localized SPs is independent of the 
size of the nanostructure Q , and that the field enhance- 
ment created in the gap between two metallic nanostruc- 
tures diverges for vanishing gap size [f|. These predic- 
tions are however in conflict both with earlier and 
with more recent experimental results, which have shown 
a size dependency of the localized SP resonance in noble 
metal nanoparticles in the size range of 1-10 nm llOl and 
pronounced deviations for dimer geometries llj, Il2j . 

This dependence of the SP resonance on the size of 
noble metal nanostructures is believed to be a signature 
of quantum properties of the free-electron gas. With de- 
creasing sizes of the nanoparticles, the quantum wave na- 
ture of the electrons is theoretically expected to manifest 
itself in the optical response due to the effects of quan- 
tum confinement [13l4l7| , quantum tunneling |l7H20| , as 
well as nonlocal response [2lM27j . Nonlocal effects arc a 
direct consequence of the inhomogencity of the electron 



gas, which arises due to the quantum wave nature and 
the many-body properties of the electron gas. 

The recent developments in analytical scanning trans- 
mission electron microscopes (STEM) equipped with a 
monochromator and electron energy-loss spectroscopy 
(EELS) [2^| gives the possibility of accessing the near- 
field energy distribution of the plasmon resonance of indi- 
vidual nanoparticles on a subnanometer scale with an en- 
ergy resolution better than 0.2 eV. This method has been 
used for the imaging of surface plasmons in many dif- 
ferent metallic nanostructures [13, [H-d^l ■ With STEM 
EELS it is possible to correlate the structural and chemi- 
cal information on the nanometer scale, such as the shape 
and the presence of organic ligands, with the spectral in- 
formation of the SP resonance of single isolated nanopar- 
ticles. STEM EELS is thus perfectly suited to probe and 
access plasmonic nanostructures and SP resonances at 
length scales where quantum mechanics is anticipated to 
become important. 

In this paper we report the experimental study of the 
SP resonance of chemically grown single Ag nanoparticles 
dispersed on 10 nm thick SisN4 membranes with STEM 
EELS. Our measurements present a significant blueshift 
of the SP resonance energy from 3.2 to 3.7 eV for par- 
ticle diameters ranging from 26 down to 3.5 nm. Our 
results also confirm very recent experiments made with 
Ag nanoparticles on different substrates using different 
STEM operating conditions [l(J, thereby strengthening 
the interpretation that the blueshift is predominantly as- 
sociated with the tight confinement of the plasma and 
the intrinsic quantum properties of the electron gas itself 
rather than ha ving an extrinsic cause. However, in con- 
trast with Ref . [lCf , we conclude that a further promotor 
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of effects due to nonlocality and quantum confinement is 
the symmetry-breaking presence of the substrate. 

We compare our experimental data to three different 
models: a purely classical local-response Drude model 
which assumes a constant electron density profile in the 
metal nanoparticle, a semiclassical local-response Drude 
model where the electron density is determined from the 
quantum mechanical problem of electrons moving in an 
infinite spherical potential well fl6| . and finally, a semi- 
classical model based on the hydrodynamic description 
of the motion of the electron gas which takes into ac- 
count nonlocal response through the internal quantum 
kinetics of the electron gas in the Thomas-Fermi (TF) 
approximation 3^, 34 1. We find good qualitative agree- 
ment between our experimental data and the two semi- 
classical models, thus supporting the anticipated nonlo- 
cal nature of SPs of Ag nanoparticles in the f-f nm size 
regime. The experimentally observed blucshift is how- 
ever significantly larger than the predictions by the two 
semiclassical models, and we attribute this difference to 
symmetry-breaking substrate effects. 



MATERIALS AND METHODS 

The nanoparticles are gro wn chemically following the 
method described in Ref. [35| and subsequently stabilized 
in an aqueous solution with borohydride ions. The mean 
size of the nanoparticles is f 2 nm with a very broad size 
distribution ranging from 3 to 30 nm. The nanoparticle 
solution is dispersed on a 10 nm thick commercially avail- 
able SiaN4 membrane (TEMwindows.com), which has a 
refractive index of approximately n ~ 2.1 [361 ]. To char- 
acterize our nanoparticles we have used an aberration- 
corrected STEM FEI Titan operated at 120 kV with a 
probe diameter of approximately 0.5 nm, and conver- 
gence and collection angles of 15 mrads and 17 mrads, 
respectively. The Titan is equipped with a monochro- 
mator allowing us to perform EELS with a energy res- 
olution of 0.15 ± 0.05 eV. We systematically performed 
EELS measurements at the surface and in the middle of 
each nanoparticle. The EELS spectra were taken with an 
exposure time of 90 ms to avoid beam damage as much 
as possible. To improve the signal-to-noise ratio we ac- 
cumulated ten to fifteen spectra for each measurement 
point. We observed no evidence of damage after each 
measurement. 

The experimental data were analyzed with the aid 
of commercially available software (Digital Micrograph) 
and three different methods were used to reconstruct and 
remove the zero-loss peak (ZLP): the first method is the 
reflected tail (RT) method, where the negative-energy 
half part of the ZLP is reflected about the zero-energy 
axis to approximate the ZLP at positive energies, while 
the second method is based on fitting the ZLP to the 
sum of a Gaussian and a Lorcntzian functions. The third 



method is to pre- record the ZLP prior to each set of EELS 
measurements. All three methods yielded consistent re- 
sults. 

The energies of the SP resonance peaks were deter- 
mined by using a nonlinear least-squares fit of our data 
to Gaussian functions. The error in the resonance energy 
is given by the 95 % confidence interval for the estimate of 
the position of the center of the Gaussian peak. Nanopar- 
ticle diameters were determined by calculating the area of 
the imaged particle and assigning to the area an effective 
diameter by assuming a perfect circular shape. The er- 
ror bars in the size therefore correspond to the deviation 
from the assumption of a circular shape, which is esti- 
mated as the difference between the largest and smallest 
diameter of the particle. 



THEORY 

In the following theoretical analysis our hypothesis is 
that the blueshift of the SP resonance energy is related 
to the properties of the electron density profile in the 
metal nanoparticle. Therefore, we use three different ap- 
proaches to model the electron density of the Ag nanopar- 
ticle. In all three approaches, we calculate the optical 
response and thereby also the resonance energies of the 
nanoparticle through the quasistatic polarizability a of a 
sphere embedded in a homogeneous background dielec- 
tric with permittivity £b- With this approach, we make 
two implicit assumptions: the first is that we can neglect 
retardation effects and the second is that we can neglect 
the symmetry-breaking effect of the substrate. We have 
validated the quasistatic approach by comparing to fully 
retarded calculations [37|], which shows excellent agree- 
ment in the particle size range we consider. The effect 
of the substrate will be taken into account indirectly by 
determining an effective homogeneous background per- 
mittivity cb using the average resonance frequency of the 
largest particles (2R > 20 nm) as the classical limit. 

The first, and simplest, approach is to assume a con- 
stant free-electron density no in the metal particle, which 
drops abruptly to zero outside the particle. This assump- 
tion produces the classical local-response Drude model 
for the response of the Ag nanoparticle, where the polar- 
izability is given by the Clausius-Mossotti relation, which 
is well-known to be size independent for subwavelength 
particles. The classical local-response polarizability 
is | 
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(1) 



where R is the radius of the particle and eo(w) = eoo(w) — 
ujp/(uj 2 + ijuj) is the classical Drude permittivity taking 
additional frequency-dependent polarization effects into 
account through eoo(w) (such as interband transitions), 



not included in the plasma response of the free-electron 
gas itself. 

The second approach is to correct the standard approx- 
imation in local-response theory of a homogeneous elec- 
tron density profile by using insight from the quantum 
wave nature of electrons to model the electron density 
profile and take into account the quantum confinement 
of the electrons. For nanometer-sized spheres, the classi- 
cal polarizability given by the Clausius-Mossotti relation 
must be altered to take into account an inhomogeneous 
electron density. In Ref. [l6| , it is shown that in general 
the local-response polarizability for a sphere embedded 
in a homogeneous material is given as 

a LQ cM = 12tt / r 2 dr / \ j_ o ' 2a 
Jo e(r,w) + 2e B 

now with a spatially varying Drude permittivity 17 1 



e(r,u) = £00(0;) - 



(2b) 



Here, n(r) is the electron density in the metal nanopar- 
ticle. Clearly, if n(r) = uq we arrive at the classical 
Clausius-Mossotti relation Eq. (fTJ) as expected. To de- 
termine the density profile in this local-response model, 
we follow the approach of Ref. and assume that the 
free electrons move in an infinite spherical potential well. 
The approach just outlined of a local-response theory 
with an inhomogeneous electron density is very similar 
to the theoretical model used in Ref. [10( for explaining 
their experimental results. It should be noted that any 
effects due to electron spill-out and quantum tunneling 
are neglected in all of the approaches that we consider. 

The third and final approach is to compare our exper- 
imental data with a linearized nonlocal hydrodynamic 
model in which the electron density is allowed to deviate 
slightly from the constant electron density used in classi- 
cal local- response theories [22|, I38l - l40| . The dynamics of 
the electron gas is governed by the semiclassical hydro- 
dynamic equation of motion |25l \26L |34| , which results in 



an inhomogeneous electron density profile. The nonlocal 
hydrodynamic polarizability onl(w) is exactly given as 



(*nl(w) = 4ttR 



epjcj) - e B (1 + £nl) 
e D (w) + 2e B (l + <W : 
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(3a) 



(3b) 



and these results constitute our nonlocal-response gen- 
eralization of the Clau sius-Mossotti relat ion of classical 
optics. Here, k^ = yuj 2 + — oS^/ e^j /3 is the wave 

vector of the additional longitudinal wave allowed to be 
excited in the hydrodynamic nonlocal theory [25ll34j . and 
ji is the spherical Bessel function of first order. Finally, 
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FIG. 1. Aberration-corrected STEM images of Ag nanopar- 
ticles with diameters (a) 15.5 nm, (b) 10 nm, and (c) 
5.5 nm, and normalized raw EELS spectra of similar-sized Ag 
nanoparticles (d-f). The EELS measurements are acquired by 
directing the electron beam to the surface of the particle. 



within TF theory (3 2 = 3/5 w|, where vp is the Fermi 



velocity [34[. We emphasize that for /3 — )• 0, the local- 
response Drude result is retrieved, since (5nl and 
Eq. (|3"a|) simplifies to the classical Clausius-Mossotti re- 
lation Eq. (TT]). 

The SP resonance energy follows theoretically from the 
Frohlich condition, i.e. we must consider the poles of 
Eq. pap . For sufficiently small blueshifts and neglecting 
damping, the resonance frequency can be approximated 
by 
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2e B P 

^/RcteooM] +2e B ' V Re M^)l 2R 



(4) 

where the first term is the common size-independent 
local-response Drude result for the SP resonance that also 
follows from Eq. (TTJ, and the second term gives the size- 
dependent blucshift due to nonlocal corrections. At this 
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stage, we note that a 1/(2R) dependence was experimen- 
tally observed in Rcfs. [y, |7| using optical spectroscopy. 
However, Eq. reveals, besides a 1/(2R) dependence, 
that there is a delicate interplay in the blueshift between 
the material parameters of the metal, through £oo(w) and 
/3, and the background medium ee- Furthermore, Eq. ^ 
shows that the blueshift can be enhanced with a large- 
permittivity background medium. 



RESULTS 

Figures QJa-c) display STEM images of Ag nanoparti- 
cles with diameters of 15.5, 10.0, and 5.5 nm respectively. 
The images show that no chemical residue is left from 
the synthesis. Thickness measurement using EELS and 
recording images of the particles at different tilt angles 
both confirmed that the overall shape is spherical, yet 
faceted (not shown). Figs.QJd-f) display raw normalized 
EELS data, acquired on similar-sized Ag nanoparticles 
with diameters 14.1, 9.8, and 6.6 nm, respectively. The 
peak corresponds to the excitation of the SP. When the 
diameter of the nanoparticlc decreases, the SP resonance 
clearly shifts progressively to higher energies. 

Figure [5] displays the resonance energy of the SP as a 
function of the diameter of the nanoparticles. A signifi- 
cant blueshift of the SP resonance of 0.5 eV is observed 
when the nanoparticle diameter decreases from 26 to 
3.5 nm. This trend is in good agreement with the results 
shown in Ref. [l(J , despite the difference in the substrate 
and the STEM operating conditions, a strong indication 
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FIG. 2. Nanoparticle SP resonance energy as a function of the 
particle diameter. The dots are EELS measurements taken at 
the surface of the particle and analyzed using the RT method, 
and the lines are theoretical predictions. We use parameters 
from Ref. u p = 8.282 eV, 7 = 0.048 eV, n = 5.9 x 

10 29 m -3 and v-p = 1.39 x 10 6 m/s. From the average large- 
particle (2R > 20 nm) resonances we determine eb = 1.53. 



that the blueshift of Ag nanoparticles is robust to ex- 
trinsic variations. Along with the EELS measurements, 
we show Eq. (Q]) for the local-response Drude model (red 
line) and the semiclassical local-response model Eq. @ 
(blue line) . Furthermore, the nonlocal relation of Eq. ^ 
(green solid line) and the approximate nonlocal relation 
of Eq. ((4]) (green dashed line) are also depicted, and we 
see that Eq. (0} is accurate for particle sizes 2R > 10 nm. 

Due to the narrow energy range in consideration 
(~ 3.0—3.9 eV), we approximate Coo(w) as a second-order 
Taylor polynomial based on the frequency-dependent val- 
ues given for Ag in Ref. [Hj]. Wc find e x (u) = (59.8 + 
i55.1)(a;/w p ) 2 -(40.3 + i42.4)(a;/wp) + (10.5 + t8.6). Since 
the refractive index of the SiaN4 substrate is hardly vary- 
ing (n w 2.1) in the narrow energy range we consider [36(, 
we assume that the background permittivity £b is con- 
stant and determine it by approximating the average res- 
onance energy of the largest particles (2R > 20 nm) as 
the classical limit, i.e. the first term of Eq. 

It is known that local Drude theory produces size- 
independent resonance frequencies of subwavclcngth par- 
ticles, but this theory is clearly inadequate to describe 
the measurements of Fig. [2] The nonlocal quasistatic hy- 
drodynamic model predicts a blueshift in agreement with 
the experimental EELS measurements. Interestingly, the 
measured blueshift is even larger than predicted. We also 
see that the local-response model with an inhomogeneous 
electron density profile shows a similar trend as the non- 
local hydrodynamic model, indicating that these two dif- 
ferent models describe very similar physical effects. The 
oscillations in the resonance energy in the inhomogeneous 
local-response model seen for small particle diameter are 
due to small variations in the density profile with decreas- 
ing size (i.e. discrete changes in the number of electrons), 
as also seen in Ref. fTp| . 

The inhomogeneous local-response model and the non- 
local hydrodynamic model, when applied to a sphere 
in a homogeneous background medium, agree qualita- 
tively with the EELS measurements, however, they do 
not provide the full picture. The issue arising is that the 
substrate is taken into account indirectly through a ho- 
mogeneous background medium, a state-of-the-art pro- 
cedure [Io| which however may not be adequate to de- 
scribe the effects of the presence of a dielectric substrate. 
It has been shown that the dielectric substrate modifies 
the absorption spectrum of an isolated sphere [42l and 
also the waveguiding properties of nanowires 
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The symmetry breaking induced by the presence of the 
substrate leads to a higher degree of electron density in- 
homogeneity in the nanoparticle, producing a stronger 
blueshift than for a particle in a homogeneous back- 
ground. The impact of the substrate on the electron den- 
sity inhomogeneity and thereby the SP resonance energy 
depends on the thickness and refractive index of the sub- 
strate, which may explain the quantitative agreement be- 
tween theory and experiment reported in Ref. (Io| . since 
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thinner substrates with smaller refractive indexes were 
used in their experiments. 

CONCLUSIONS 

We have investigated the surface plasmon resonance of 
spherical silver nanoparticlcs ranging from 26 down to 
3.5 nm in size with STEM EELS and observed a signif- 
icant blueshift of 0.5 eV of the resonance energy. We 
have compared our experimental data with three differ- 
ent models based on the quasistatic optical polarizability 
of a sphere embedded in a homogeneous material. Two 
of the models, a nonlocal hydrodynamic model and a 
generalized local model, incorporate the inhomogeneity 
of the electron density induced by the quantum wave na- 
ture of the electrons. These two different models produce 
similar results in the SP resonance energy and describe 
qualitatively the blueshift observed in our measurements. 
Although our exact hydrodynamic generalization of the 
Clausius-Mossotti relation predicts a nonlocal blueshift 
that grows fast [as 1/(2R)} when decreasing the diam- 
eter and increases even faster for the smallest particles 
(2R < 10 nm), the observed blueshifts are nevertheless 
larger than predicted. 

The quantitative agreement between the two different 
theoretical models and the discrepancy with the larger 
observed blueshift suggest that the assumption of a ho- 
mogeneous background is inadequate to describe the 
presence of the substrate. In contrast to the quantitative 
agreement between theory and experiment reported in 
Ref. we conclude that state-of-the-art models with 
a homogeneous dielectric background permittivity can- 
not explain the sizable resonance blueshifts that we ob- 
served in our EELS measurements. Further theoretical 
and experimental investigations are needed to quantita- 
tively understand the effect of the substrate on the elec- 
tron density inhomogeneity for metallic nanoparticlcs be- 
low ten nanometer in size. 
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